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Tully-Fisher relation, key to dark companion of baryonic matter 

Yousef Sobouti, * Akram Hasani Zonoozi, and Hosein Haghi *** 

Institute for Advanced Studies in Basic Sciences (lASBS), P. O. Box 45195-1159, Zanjan, Iran 
Received Accepted ***** 

Abstract. Rotation curves of spiral galaxies i) fall off much less steeply than the Keplerian curves do, and ii) 
have asymptotic speeds almost proportional to the fourth root of the mass of the galaxy, the Tully-Fisher relation. 
These features alone are sufficient for assigning a dark companion to the galaxy in an unambiguous way. In regions 
outside a spherical system, we design a spherically symmetric spacetime to accommodate the peculiarities just 
mentioned. Gravitation emerges in excess of what the observable matter can produce. We attribute the excess 
gravitation to a hypothetical, dark, perfect fluid companion to the galaxy and resort to the Tully-Fisher relation 
to deduce its density and pressure. The dark density turns out to be proportional to the square root of the mass 
of the galaxy and to fall off as r~''^+™', a <C 1. The dark equation of state is barrotropic. For the interior of the 
configuration, we require the continuity of the total force field at the boundary of the system. This enables us to 
determine the size and the distribution of the interior dark density and pressure in terms of the structure of the 
observable matter. The formalism is nonlocal and nonlinear, and the density and pressure of the dark matter at 
any spacetime point turn out to depend on certain integrals of the baryonic matter over all or parts of the system 
in a nonlinear manner. 
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1. Introduction 

Gravitation of the observable matter in galaxies and clus- 
ters of galaxies is not sufficient for explaining their dynam- 
ics. Dark matter scenarios and/or alternative theories of 
gravitation (see e. g., Milgrom 1983; Behar and Carmelli 
2000; Capozziello et al. 2002, 2003 and 2006; Carroll et al. 
2004; Norjiri et al. 2003 and 2004; Moffat 2005; Sobouti 
2007) are called in to resolve the dilemma. The fact re- 
mains, however, that the proponents of dark matter have 
always looked for it in observable matter. No one has, 
so far, reported a case where there is still no baryonic 
matter, but there is a dynamical issue to be settled. In 
view of this negative observation, it has been conjectured 
( [Sobouti 2008 a, b; 2009[ ) that, if there is a dark compan- 
ion to any baryonic matter, there must be rules to connect 
the properties of the twin entities. On the other hand, the 
existence of such a rule will entitle one to interpret the 
case as an alternative gravity, thus reducing the difference 
between the two paradigms to the level of semantics. This 
conclusion, however, is true as long as the assumed dark 
matter does not interact with the baryonic one in any 
other way than through its gravitation. 
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Sobouti assumes a spherically symmetric system, at- 
tributes a dark perfect fluid companion to it, and re- 
quires the rotation curve of the system to display the same 
asymptotic behavior as those of the actual spirals. The 
reason for the assumption of a dark fluid instead of the 
conventionally assumed dark pressureless dust, is to en- 
sure the satisfaction of the Bianchi identities and thereby 
the baryonic conservation laws (See sect. 8 for further ex- 
planation) . 

In regions outside to the baryonic system, he finds the 
density and pressure of the dark fluid companion in terms 
of the mass of the host system. The Tully-Fisher rela- 
tion and the slow non-Keplerian decline of the rotation 
curves play key roles in determining the relation between 
the matter and its dark twin. 

In this paper, we follow the same line of argument to 
find the structure of the dark matter in the interior of the 
baryonic system. The continuity of the total gravitational 
force at the boundary of the observable matter leads to the 
dark matter distribution in the interior. The Tully-Fisher 
relation is a nonlocal and nonlinear feature of the dynam- 
ics of galaxies: a) The presence of the total or partial in- 
tegrals of the baryonic matter in the structure of both 
exterior and interior solutions reflects the nonlocality. b) 
That the excess gravitation does not increase proportion- 
ally upon increasing the mass of the host galaxy indicates 
the nonlinearity. To emphasize these two features, we re- 
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fer to the formahsm developed here as the nonlocal and 
nonlinear (NN) one. 

To check its validity, the formalism is applied to NGC 
2903 and NGC 1560, two examples of high and low sur- 
face brightness galaxies, respectively, and the resulting ro- 
tation curves are compared with those obtained through 
other approaches. 

2. Model and formalism 

The following is a brief background from Sobouti (2008a, 
b; 2009). The physical system is a spherically symmetric 
baryonic matter of finite extent. By conjecture there is a 
dark presence that pervades both the interior and exterior 
of the system. The spacetime metric inside and outside of 
the system is necessarily spherically symmetric and takes 
the form 



ds^ = ^B{r)dt^ + A{r)dr'^ + r^{de^ + sin^ 9dip^). 



(1) 



Let both the baryonic matter and its dark companion be 
perfect fluids of densities p, pd, of pressures p, pd, respec- 
tively, and be at rest. From the field equations of general 
relativity (GR), we find 



1 



d_ 

dr 



- 1 



= -(P + Pd), 



+ ^ = Hp + Pd) + {p + Pd)], 



(2) 



(3) 



where we have let 8ttG — = 1, and "'— d/dr. In the 
nonrelativistic regime, we neglect the pressures, eliminate 
the densities between the two equations, and arrive at 



B' I, . 

-^ = -{A-l). 

B r 



(4) 



In the following two sections we solve Eqs. (H]) - (21) inside 
and outside the baryonic system. 

3. Exterior solution 

Hereafter, the parameters pertaining to the interior and 
exterior of the system will be labeled by the superscripts 
(i) and (e), respectively. The unknowns in Eqs. are 
A, B, Pd, Pd, and the dark equation of state. We begin 
with Eq. Q and assume that in the baryonic vacuum, 
p = p = 0, the factor (A^*^' — 1) is differentiable and has 
the series expansion 



r>R, 



where the indicial exponent a and sq are dimensionless, 
si has the dimension of length, tq is an arbitrary length 
scale of the system, and R is the radius of the baryonic 
sphere. Substituting Eq. ^ into Eq. (HI and integrating 
the resulting expression, gives 



exp 



So si 
a (1 + a)r 



(6) 



We expand the exponential factor, keep its first two terms, 
and for the weak field gravitational potential, (j) = {B — 
l)/2, find 



a 



(1 -f- a)r 



(7) 



The square of the circular speed of a test object orbiting 
the galaxy is 



dr 



.0 + ^ 

r 



(8) 



Equation ^ is the rotation curve of our hypothetical 
galaxy in its baryonic vacuum. It has an asymptotically 
constant logarithmic slope, 

A = dlnv^ /dlnr —a as r ^ cxd. 



3.1. Determination of a, sq, si, ••• 

Rotation curves of actual spiral galaxies have two distinct 
non-classical features: 

i) Their asymptotic slopes are much flatter than 
that of the Keplerian curves, —1, (|Sanders 1996 



Bosma 19811 |Be gmann 1989[ Persic and Salucci, 1995| 



Begmann et aL 1991] |Sanders and Verheijen 19981 



Sanders and McGhaugh 2002). This implies a <C 1. From 
Persic et al. 1996i who study 1100 galaxies with the aim 
of arriving at a universal rotation curve, we estimate 



a < 0.01. 



(9) 



Moreover, a does not seem to be a universal constant. The 
rotation curves of more massive galaxies appear to fall off 
somewhat more steeply than those of the less massive ones 
(jPersic et al. 1996| . Hereafter, for simplicity but mainly 
for pedagogical reasons, we work in the limit of a 0. 

ii) Their asymptotic speeds follow the TuUy-Fisher re- 
lation. They are almost proportional to the fourth root 
of the mass of the host galaxy ( TuUy and Fisher 1977[ 
IBegmann 1989} |McGaugh et al. 2000| |McGhaugh 2005| ). 
In Eq. ([SI), letting a — > 0, the dominant term at large 
distances is = so/2. We identify this v with the TuUy- 



Fisher asymptote and conclude that 
So = A [M/Mq)^''^ , A = 2.8 X IQ-^' 



(10) 



where M is the galactic mass, and A can be obtained ei- 
ther from a direct examination of the observed asymp- 
totic speeds (jSobouti 2007|) or from a comparison of the 
first term of Eq. ^ with the low acceleration limit of 
MONO dMilgrom 1983|: v^/r {aogn)^'^, oo = 1.2 x 



(5) 10-i°m sec-2 pegmann 1989] ) 



Again letting a — ^ 0, the second term in Eq. ([8]) is 
the classic Newtonian or GR term. Therefore, si should 
be identified with the Schwarzschild radius of the host 
galaxy: 



si = 2GM/c^ 



(11) 



Here, for clarity, we have restored the constants and G 
and written si in physical units. There is no compelling 
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observational evidence to indicate the need for other terms 
in Eqs. ([5]) - ([8]). Therefore, at least at the present state 
of the extent and accuracy of the observational data, we 
truncate the series at the si term. 

4. Interior solution 

The first and foremost condition to be satisfied is the con- 
tinuity of the total force exerted on a test object at the 
boundary, R, of the baryonic system. Pressure forces are 
anticipated to be insignificant in the present problem so 
are ignored. The gravitational forces remain. From Eqs. 
© - ((TT|) . the exterior force is 



dr 



A 



M 

Mo 



1/2 



f 2GM 1 

^ 2 9 

r c r 



r> R. 



(12) 



By analogy, for the interior of the system we adopt 



dr 



A 



2GM{r) 1 



, r < i?, (13) 



where M{r) — At: pr^dr is the variable baryonic mass 
inside the radius r. The continuity of the exterior and 
interior forces at the boundary is evident, QED. Once 
the baryonic p{r) and Af(r) are known, 0(')(r), B^{r) « 
1 + 2<^(^) and can be integrated. The expression for 
the latter is much simpler and is given below for later 
reference. From Eq. (jj]) we find 



dr 



( M{r 



1/2 



2GM{r) 1 



(14) 



This has the same form as Eq. ([S]), where M is replaced 
by M{r). 

5. Structure of the dark matter 

The densities are obtained from Eq. ^ or equivalently 
from Poisson's equation through Eqs. - For the 
exterior dark density we find 



X 



M 



1/2 



r>R. 



(15) 



Note the square root dependence on the mass of the galaxy 
and the fall out as r~^. For the interior, A'*^ is given by 
Eq. (I14|) . whose first term gives the interior dark density 
and the second renders the baryonic density, p. Thus, 



(r)=A 



"Af(r)' 


1/2 ^ 







1 + 27r 



M{r) 



The dark matter inside the radius r is 

"M(r) 



Mrf(r) = An 



i(r)r dr = X 



Mr.- 



1/2 



r<R. (16) 



(17) 



Equation (flT)) holds for any r. For r > R, however, M(r) 
attains its maximum constant value, M, and Mij{r > R) 
becomes proportional to r. 



It is instructive to look at the behavior of Eq. pB|) in 
the neighborhood of the origin, where p — > pc and M{r) — > 
Airpcr^/S. Equation (fTH)) tends toward 



n^ ^ \ f ^'^ P'' \^'^ -1/2 
5y(r 0) -A I — I r / 



2 V 3 A/0 



Similarly, 



Mf(r-.0) = A('i^^ 
^ ' I 3 Af0 



1/2 



^5/2 



(18) 



(19) 



While the density becomes singular as r 0, no cusp 
develops. For the measure r'^dr tends to zero as r — > 0. 

Pressures of the matter and of its dark companion 
are obtained from their hydrostatic equilibrium, a require- 
ment of the Bianchi identities. The general formula is 



V_ 

P 



2 S" 



dr 



(20) 



For the exterior pressure from Eqs. ((20|) . (fT5|) . (fT2|) . we 
find 



p^;\r) 



1 



fo 



2 si 
3^ 



r> R. 



(21) 



The presence of an extra factor of sq in Eq. (j2T|l makes 
the pressure an order of magnitude less than the density 
and justifies the approximation made in the derivation 
of Eq.(|4]) and thereafter. The equation of state, p{p), in 
the exterior region is obtained by eliminating r between 
Eqs. ((2T|) and (fT5| . It is barrotropic. The internal pressure 
is obtained in a similar way. It is, however, an involved 
expression so is too involved expression to give here. 

A pedagogical note: Throughout the text, except in 
Eg. pTj) . we have chosen SttG = = 1. To write the 
results in physical units, the rule is to multiply, every- 
where, the potentials, (j), by c^, the dark densities, pd, dark 
masses, Af^, by c^/SttG, and the dark pressures, pd, by 
cVSttG. 

6. Application to actual spirals 

Spiral galaxies are flattened objects. Their approximation 
as spherical systems introduces an error on the order of 
[Rgyr/rY, where Rgyr{r) is the gyration radius of the 
mass enclosed within a radius r. In a flat system that 
thins out as an exponential or as a Kuzmin disk, say, this 
ratio would be a few parts in thousand and small enough 
for our purpose. This is also the practice of all the authors 
quoted so far in this paper. To illustrate the practical ap- 
plicability of the formalism developed here, we construct 
the rotation curves of two standard high- and low- surface 
brightness galaxies and compare the results with those 
obtained through MOND's formalism. 

NGC 2903 is a textbook example of a high surface 
brightness spiral. It has a large stellar component and 
small HI content. The gas is confined to the galactic plane 
and follows circular orbits. It is well observed out to about 
40 kpc ( iBegmannT989l ). In contrast, NGC 1560 is a low 
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Table 1. Minimum ^ and fitted stellar mass-to- light ra- 
tio, T, of MOND and of our NN formalism. 



Galaxy 



X^iV TiVJV xliOND Tmond 



NGC 2903 
NGC 1560 



4.97 
1.52 



1.7 
0.3 



6.07 
3.35 



3.0 
1.1 



surface brightness spiral with a dominant gas component. 
Its observed rotation curve extends out to about 8 kpc 
and does not seem to have reached its asymptotic regime. 

In Fig. 1 we construct the rotation curves of our NN 
formalism from Eq. ([T^ . in which M(r) is the total, stellar 
plus HI, mass interior to r. The free adjustable parameter 
in matching the theoretical curves to data points, is the 
'stellar' mass-to-light ratio, T, assumed to be constant 
throughout the galaxy. For comparison we have also in- 
cluded the rotation curves of MOND. That the NN curves 
trace the data points more closely than the MOND ones 
can be seen pictorially. The test and T's of Table 
1, however, illustrate this in a quantitative way. In both 
galaxies our x'nn noticeably small. Significant, however, 
is the low stellar mass-to-light ratio of the young and gas- 
dominated NGC 1560. Our T — 0.3 is, by far, closer to 
0.4 estimate of McGaugh (2002) than to 1.1 of MOND. 

Our next project is to study pressure-supported sys- 
tems, globular clusters and dwarf spheroidal galaxies 
(dSph). Globular clusters are commonly believed to be 
almost Newtonian systems, while dSph's show significant 
deviations from Newtonian regimes. The low baryonic 
mass and extremely high dynamical-mass-to-light ratio of 
dSph's are inconsistent with population synthesis mod- 
els (jHilker 20061 IJordi et al. 20091 [Angus 2008] ). Our ap- 
proach is to find a counterpart of the classical virial the- 
orem for our proposed gravity and to solve a modified 
Jeans equation. The aim is to verify whether the velocity 
dispersions obtained via Jeans equation fit the observed 
data. We also hope to be able to come up with a notion 
equivalent to the fundamental plane for galaxies where 
one arranges the galaxies on a two-parameter-plane in a 
three-dimensional space of luminosity, velocity dispersion, 
and some other global characteristics of the galaxies. 
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Fig. 1. Points with error bars are observed data. Dotted 
and dashed lines are the contributions of the gaseous 
and stellar components to the rotation curves, respec- 
tively. Dashed-dottcd line is the rotation curve con- 
structed through MOND's formalism. Solid line is our ro- 
tation curve calculated from Eq. (fT3)) . The free parameter 
in matching theoretical curves to data points, is the stellar 
mass-to-light ratio. 



7. Nonlocality and nonlinearity of the formalism: 

The masses M and A/(r) are integrals over all or parts of 
the system. Their presence, in the structure of the space- 
time metric, in the rotation curve, in the expressions for 
the dark densities and pressures, etc., reflects the nonlo- 
cal nature of the theory. That these integrals enter the 
formalism not in a linear way indicates the nonlinearity of 
it. Both features are rooted in the Tully-Fisher relation, 
which requires the dynamical variables at one spacetime 
point to depend on the integral properties of the whole or 
parts of the system through the square root of these inte- 
grals. Any attempt to derive the spacetime metric enter- 



tained in this paper through a variational principle should 
take these two features into account. 

In this respect, Hehl and Mashhoon's generalization 
of GR, ( |Hehl and Mashhoon 2009 a, b[ ), constructed 
within the framework of the translational gauge theory 
of gravity, is interesting. In the weak field approximation, 
the excess gravitation coming from the nonlocality of 
their theory can be interpreted as a dark companion to 
the baryonic matter. In the case of a point baryonic mass, 
M, the dark density has the expected distribution, 
But it does not obey the TuUy - Fisher relation. Instead 
of M^/^, it is proportional to Al itself. 
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8. Concluding remarks 

The formahsm developed here is a dark matter scenario 
or, equivalently, a modified GR paradigm to understand 
the non-classical behavior of the rotation curves of spi- 
ral galaxies. Following ( [Sobouti 2008 a, b; 2009P , we at- 
tribute a hypothetical dark perfect fluid companion to our 
model galaxy , and find the size and the distribution of 
the companion by comparing the rotation curve of the 
model with those of the actual galaxies. However, as long 
as the dark companion displays no physical characteristics 
other than its gravitation, one has the option to interpret 
the scenario as an alternative theory of gravitation. Here, 
for example, one may maintain that the gravitation of a 
baryonic sphere is not what Newton or Schwarzschild pro- 
fess, but rather what one infers from the spacetime metric 
detailed above. In fact we wish to emphasize that any 
modified gravity is expressible in terms of a dark matter 
scenario. And vice versa, any dark matter paradigm, in 
which the matter and its dark twin are related by certain 
rules, is explainable by a modified gravity. The difference 
between the two alternatives is semantic. 

Dynamics of galaxies is a nonrelativistic issue. Yet, its 
analysis in a GR context answers questions that otherwise 
are left out. In particular, in a nonrelativistic scenario, 
there is neither need nor logic to assign a pressure field 
to a hypothetical matter that one knows nothing about 
its nature. In a GR context, on the other hand, the dark 
matter has to have a pressure field and has to be in hy- 
drodynamic equilibrium as a requirement of the Bianchi 
identities and thereby of the conservation laws of the bary- 
onic matter, i.e. the vanishing of the 4-divergence of both 
sides of the field equations. Let us also note in passing 
that all those metric approaches that attempt to explain 
the galaxy problems with the aid of a single scalar field 
are subject to the same criticism, namely the violation of 
the Bianchi identities and of the conservation laws. 

Regions outside to the baryonic matter are not dark 
matter vacua. Therefore, the Ricci scalar does not vanish, 
and there are excess lensing and excess periastron preces- 
sion caused by the dark matter. These are analyzed in 
Sobouti (2008a, b; 2009). 

The formalism is good for spherical distributions of 
baryonic matters. An axiomatic generalization to non- 
spherical configurations or to many body systems re- 
quires further deliberations and more accurate observa- 
tional data to help find some solutions. One might need 
other postulates not contemplated. The difficulty lies in 
the nonlinearity of the formalism. There is no superpo- 
sition principle. One may not add the fields of the dark 
companions of two separate baryonic systems because sq 
of Eq. ([ini) is not fincar in M or M(r). 
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